Load data

Load data

Data preparation

Extract Harvest & plant dates into months and years

##      cname  ccode   planting  lai gwad    harvest h_month h_year p_month p_year
## 1 Alebtong UG3067 1990-04-22 0.15 1368 1990-09-20       9   1990       4   1990
## 2 Alebtong UG3067 1990-04-22 0.15 1369 1990-09-20       9   1990       4   1990
## 3 Alebtong UG3067 1990-04-23 0.22  542 1990-09-20       9   1990       4   1990
## 4 Alebtong UG3067 1990-04-23 0.15 1266 1990-09-02       9   1990       4   1990
## 5 Alebtong UG3067 1990-04-23 0.16 1356 1990-09-20       9   1990       4   1990

Let’s define/extract the various seasons in Kenya.

lets now aggregate RHEAS estimated yields per year per county.

Plots

Annual maize yields variation.

Country level yields.

Display maps

Load country shapefile

## Loading required package: sp
## Warning: multiple methods tables found for 'approxNA'
## tmap mode set to interactive viewing
## Warning in sp::proj4string(obj): CRS object has comment, which is lost in output
## Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1

Merge boundary with average distict level yield information.

Visualize average district level information.

## tmap mode set to interactive viewing
## Warning in sp::proj4string(obj): CRS object has comment, which is lost in output